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Breaking time-reversal symmetry (TRS) in the absence of a net bias can give rise to directed 
steady-state non-equilibrium transport phenomena such as ratchet effects. Here we present, the¬ 
oretically and experimentally, the concept of a Lissajous rocking ratchet as an instrument based 
on breaking TRS. Our system is a semiconductor quantum dot (QD) with periodically modulated 
dot-lead tunnel barriers. Broken TRS gives rise to single electron tunneling current. Its direction 
is fully controlled by exploring frequency and phase relations between the two barrier modulations. 

The concept of Lissajous ratchets can be realized in a large variety of different systems, including 
nano-electrical, nano-electromechanical or superconducting circuits. It promises applications based 
on a detailed on-chip comparison of radio-frequency signals. 

PACS numbers: 73.63.Kv, 05.60.-k, 85.35.Gv 


Ratchets cause directed particle motion due to a com¬ 
bination of broken symmetry and non-equilibrium forces, 
where the latter may be deterministic or fluctuating. 
The most famous example is Feynman’s flashing ratchet 
which uses a pulsating spatially asymmetric potential 
to actively turn fluctuations into work urn- Another 
species is the rocking ratchet driven by forces periodic in 
time but with broken spatio-temporal symmetry. A sim¬ 
ple example of a rocking ratchet is a pump which trans¬ 
ports electrons one-by-one through a QD driven by two 
external periodic forces with a relative phase breaking 
the symmetry EE]. This is in contrast to the somewhat 
simpler turnstile where the spatial symmetry is broken by 
a finite dc voltage [7] . In the non-adiabatic limit electron 
pumps are investigated for their suitability as current 
standard m- 

In this article we restrict ourselves to the adiabatic 
regime and study a generic implementation of a rock¬ 
ing ratchet by applying two time-periodic forces, phase- 
locked at various commensurate frequencies. In our im¬ 
plementation we measure the dc current I through a 
QD embedded in the two-dimensional electron system 
(2DES) 90 nm beneath the surface of an etched 1 /im wide 
channel of GaAs/AlGaAs heterostructure. The 2DES 
is cooled to ^lOOmK where its carrier density is Ue — 
2.83 X 10^^ m“^ and its mobility is jie — 320 m^ s“^. 
We control the QD by applying voltages to two metal 
gates [see lower inset in Eig. in»)i' 

Eor a first orientation we present in Eig. [^a) a sta¬ 
bility diagram of our QD measured at finite dc voltage 
V = (hR ~ = 100 juV applied between its two leads 

(at chemical potentials /iL,R)- Plotted is the current I 
as a function of gate voltages Vl and Vr applied to the 
left (L) versus right (R) gate [yellow in the rhs inset of 
Eig. [h^)]- The axes El,r are offset relative to the ac¬ 
tually applied voltages, such that W = hp = 0 at our 
working point, which is marked in Eig. [^a) by a black 


cross. We define the chemical potential of fho QD as 
the energy needed to add the next electron to it, where 
n is an index number and we choose n = 0 for the dot 
level closest to the working point. Eor /ip ^ hn ^ Pl tPo 
Goulomb blockade is lifted and a current I fiows along 
broadened lines oriented perpendicula£ to the bisecting 
line (dashed line). With increasing Ll,R 5 the barriers 
coupling the QD to both leads shrink causing I to grow. 
At the bisecting line, the QD is symmetrically coupled 
to both leads but with increasing distance, one barrier 
shrinks and the other one grows, causing I to decrease. 

We parametrize the QD by its dot-lead tunnel rates 
rL,R and its charging energies Ec{n) = /in where 

at our working point, Ec(0) = 1.27meV. Explor¬ 

ing a limited range of the stability diagram the capac¬ 
itive coupling between gates and the QD^can be mod¬ 
eled by the linear relation /i^ = /i^ — + Lp) with 

0.38e. Here e is the elementary charge and /i^ are 
the chemical potentials for Lp = Hr = 0. The tunnel cou¬ 
plings rL(hL) and rp(LR) are controlled by the respec¬ 
tive gate voltages Ll,r, where we neglect the very weak 
cross couplings rL(T4) and rp(EL)- The rL,R(EL,R) are 
determined by fitting a theoretical prediction to the cur¬ 
rent in Eig. [^a). Our calibrations are detailed in the 
Supplemental Material, Ref. [10]. In the following we 
define rocking ratchets by making use of the mutually 
coupled dependencies of rL,R and /i^ on both gate volt¬ 
ages. 

To break the TRS, we periodically modulate Vi,{t) and 
Lp(t) such that they define trajectories along Lissajous 
figures centered at the working point: 


Vk,4>{t) 


= A f cos(27r/i - (/))\ 
\VR{t)j V cos(fc27r/i) ) ’ 


( 1 ) 


where, for simplicity, we restrict ourselves to integer fre¬ 
quency ratios k. We prepare the two modulation am- 


2 


a 



Vr (mV) 




Vr (mV) 


Vr (mV) 


FIG. 1. (a) Current I through the QD as function of gate voltages 14 ,r in response to a dc voltage V = (/xr —/iL)/e = 100 juV. 

A double arrow indicates the charging energy at the working point (black cross). Lower inset: scanning electron microscope 
image of the sample. The wafer is etched at the gray area, only the blue channel contains 2DES. Two Ti/Au gates (yellow) 
are used to define a QD, a third gate (gray) is electrically grounded. Upper inset: QD sketch; vertical lines are electrostatic 
barriers with tunnel couplings rL,R controlled by gate voltages. Horizontal lines are chemical potentials fin of the QD, blue 
areas indicate occupied states in the degenerate leads at chemical potentials /xl,r. (b) Ul,r are modulated according to Eq. 0’ 
both with frequency / = 200 MHz (/c = 1). Measured (left) versus calculated (right, V = 1 A^V) current I versus modulation 
amplitude and phase. Inset: measured /(/) at A = 9mV for the two phases 0 ~ O.Stt (blue), 0 ~ O.Gtt (red), (c) Vector field 
a (explained in text and Ref. [TO]) versus gate voltages and typical Lissajous trajectories for /c = 1,2,3 according to Eq. 0 . 
Black crosses mark the working point, QD sketches indicate relative strengths of tunnel couplings to the leads. 


plitudes to be identical, A. The trajectories can cross 
several charging lines of the stability diagram. 

Our first experiment resembles an adiabatic electron 
pump defined by Eq. 0 with k = 1, i. e. identical fre¬ 
quencies. Figure [^b) compares a typical measurement 
(Ihs) with model predictions (rhs), both for / = 200 MHz 
and /ip — k'R- The current, plotted as a function of A 
and phase 0, displays clear resonances following lines of 
constant Acos(0/2) (solid lines). These resonances indi¬ 
cate discrete current contributions of the n-electron QD 
levels separated by n^/AA = Ec{n)/2, where the fac¬ 
tor 2 accounts for two gates being modulated. During 
each pumping cycle approximately mt{2 AK^/Eq) ground 
state levels fulfill = A^l,r twice. For / < 200MHz 

we find / oc / at fixed 0 and A [inset in Fig. j^b)] which 
confirms that we operate in the adiabatic regime, where 
the QD is always in its momentary ground state and the 
charge transferred across it per cycle is constant m- 
Each QD level is unoccupied whenever /in(t) > and 
its n-electron ground state is occupied for /in(t) < /iL,R- 
Whenever En{t) ^ /iL,R an electron tunnels into the QD 
if d/j^n/dt < 0 and out of the QD if d/i^/dt > 0. With 
which lead the electron is thereby exchanged depends 
on (p and the ratio Tl/Tr, which is modulated in time, 
see Eq. 0 and above. 

This scenario resembles a rocking ratchet [4] where the 
current changes direction at 0 = 0, tt, independent of k. 
For k = 1 it gives rise to the two-fold symmetry of a 
pump observed in Fig.j^b). The symmetry can be com¬ 
promised for four reasons (independent of k): (i) dissipa¬ 
tion by transitions within the QD’s excitation spectrum. 


This, however, would go along with non-adiabaticity 
which we already excluded [inset in Fig. [^b)]; (ii) a 
spatially asymmetric local disorder potential which in¬ 
fluences the ratio Tl/Tr as function of n, depending on 
the electronic probability distribution. This is ignored in 
our model but can explain deviations to the measured 
current, particularly near 0 = 0 where the QD’s chem¬ 
ical potentials are modulated strongest and addi¬ 

tional spatial disorder causes rectification El; (hi) the 
n-dependence of the charging energy and the choice of 
the working point; (iv) a dc voltage between the leads 
{V 7 ^ 0). Points (hi) and (iv) are included in our model 
and discussed below. 

A formal description of the general case for any inte¬ 
ger k including different frequencies modulating the two 
gates, is detailed in Ref. m and can be summarized 
as follows: We start from an expression for the current 
through the time-periodic system obtained with Floquet 
transport theory mun]. Taylor expansion up to first or¬ 
der in Q provides the adiabatic limit of the dc current: 

/ = GU + (2) 

which can be traced back to the transport properties of 
the time-independent system for the parameters along 
the trajectory f4^0(t), defined in Eq. 0 . The first 
term of Eq. 0 contains the time-averaged conductance 
of the QD, G. It can be understood as dc current 
which flows whenever any resonance fulfills the condition 
Ml < Mn ^ Mr (or Mr ^ Mn ^ Ml). For the theory data 
in Fig. 0;b) we have used V = 1 fiY which provides good 
agreement with the measured data. This tiny dc volt- 































3 



FIG. 2. Measured (left) and computed (right) current as a function of modulation amplitude A and phase difference 0 between 
Vh{t) and VK{t): (a) I4(t) at 50 MHz and VR{t) at 100 MHz {k = 2), see Eq. Q, and (b) at 50 MHz and 150 MHz {k = 3). 
Solid lines indicate the resonances fin — IJ^l.r for n — 3,4. Bottom: Lissajous figures for k — 2 and /c = 3. 


age is not purposely applied but caused by the current 
amplifier. 

The second term of Eq. Q is the dc current caused 
by the pumped charge, . This ratchet current can 

be expressed as a closed loop integral of a vector field 
a{v) [v being defined in Eq. Q], which follows from the 
scattering matrix m- Here, we restrict ourselves to a 
qualitative ^planation and plot its component ai as a 
function of Vl,r in Fig. j^c); note that \a 2 \ — \ai\. a{v) 
vanishes unless inL,R and, hence, marks the charg¬ 

ing lines. Its direction relates to the ratio Tl/Tr: since 
ai ±a 2 , a is either parallel or anti-parallel to the (1,1) 
direction in Fig. [^c), depending on whether Vh ^ Vr 
(implying Tp > Tr) or Ft < 

The third contribution in Eq. ([^, stems from 

a tiny ac modulation of the lead chemical potentials in¬ 
duced by the capacitive coupling between each gate and 
its adjacent lead. The amplitude of the resulting ac 
source-drain voltage is of the order 0.1/iV and, owing to 
the capacitive coupling, is phase-shifted by —7r/2 with re¬ 
spect to the gate voltages. This phase shift between the 
ac bias and the time dependent QD conductance gives 
rise to dynamic rectification. turns out essential for 
the quantitative agreement between the measured cur¬ 
rents and the theoretical predictions of our model. Inter¬ 
estingly, 7 ^ 0 requires a phase shift between the two 
modulated gate voltages. As such, dynamic rectification 
is very different in nature than the usually and hitherto 
discussed rectification caused by static asymmetries of 
the geometry and I-V curve M- 

The symmetries of these three contributions as a func¬ 
tion of (j) can be revealed by formally reverting time [ID]. 
We find 2k symmetry points at which G obeys TRS, while 
Qcycie jrect anti-symmctric under time reversal 
and, hence, under phase reflection. 

Figure l^c) contains example Lissajous figures Vk,(p{t) 
for (p = 7r/{k 1): a circle for /c = 1, a distorted fig¬ 

ure eight for k = 2 and a triple loop for k = 3. For 


k = 1 and p slightly different from tt, Vk=i,(f){t) is an 
eccentric ellipse which corresponds to the pumping mea¬ 
surements already discussed in literature [5l [6] [dashed 
line in Fig. [^c)]. Whenever Vk,(i){t) crosses a charging 
line, the charge of the QD changes by one electron. If 
this happens in a red region corresponding to ai,a 2 > 0 
with Fp > Tr an electron will be exchanged preferably 
with the left lead. In the blue region with ai,a 2 < 0 
and Fp < Fr, charge exchange with the right lead is pre¬ 
ferred. We define / > 0 for electrons fiowing from the 
right to the left. Then, crossing a blue area from above 
and a red area from below (see example for k = 1) both 
contribute to / < 0, each with half an electron charge 
per cycle. Our example for k = 3 also results in / <0. 
For /c = 2, the contributions of the outer loops to 
have opposite sign and, thus, cancel each other to some 
extent, despite that the curves have a different symme¬ 
try axis than the resonance lines. The same holds for 
the inner loops. Therefore, we expect the pump current 
for k = 2 to be generally smaller than for /c = 1. For 
/c = 3, by contrast, symmetry-related parts of each curve 
again contribute to with the same sign. A general¬ 

ization of these arguments leads to the expectation that 
for even values of /c, the pump current should be larger 
than for odd values. For larger /c, however, there will be 
an increasing number of contributions with any sign and, 
thus, the situation becomes less clear. 

In a realistic system, in addition, the symmetry of 
a{v) is compromised because of point (iii) above: the 
separation of charging lines, Eq, depends on n and our 
working point (see Fig. is chosen to further decrease 
symmetry. As a result the contributions to from 

the outer loops for k > 1 do not completely cancel each 
other and we find a ratchet current for all k, specifically 
with deviations from integer for k = 3, 5, 7,... 

and from = 0 for k = 2,4,6,.... The mapping 

(j) ^ —(j) should survive the compromised symmetry dis¬ 
cussed above and still result in a reversal of the direc- 
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FIG. 3. Symbols are measured, lines depict theory, (a) 7(0) near zero bias for ~ 27^2^, I/off — — 1/iV (blue) and 
~ 0.57^^, I4ff ~ 0/iV (red); k — f — 100 MHz. (b) 7(0) for ti^A ~ 27^^, I4ff ~ 9 /xV including curves at finite bias. 
Arrows indicate the phases at which the QD levels with indices n start to reach the transport window (/xl < /Xn < Mr or 
Mr ^ Mn < Ml) during each modulation period. The relative shift between n — 0,-1 and n — 1,-2, etc. is caused by our 
working point not being centered between two QD chemical potentials [see Fig.[^a)]. (c) 7-I/eff characteristics for ~ 27^^, 
at current peaks near black arrows in (b). The maxima slightly shift proportional to |I4ff|. Upper right inset: magnihcation, 
theory lines intersect at I4ff ~ 9 mV. Upper left inset: average conductance G. Lower inset: (red, negative values) and 

pump current (black, positive), both being independent of Ueff- 


tion of Vk,cf){t) and, therefore, also in a current reversal. 
Furthermore, the current direction should posses the k- 
fold symmetry 0^0 + 27r//c. As a pump has two-fold 
symmetry, the case of k > 1 goes qualitatively beyond 
the scope of a pump. For the 2k equally spaced phases 
0 = 0, tt/A:, ..., {2k — l)7r/7, no charge is pumped at all, 
Qcycie should change its sign and, thus, vanish. This 
symmetry is apparent albeit not perfect in Fig. which 
compares the measured ratchet currents for A: = 2 and 
A: = 3 as function of the modulation amplitude and 0 with 
model calculations. More values of k are available in Ref. 
m- To improve modeling of our measured data we in¬ 
cluded all experimentally known facts such as Eq (n) and 
line broadening into the numeric calculations (see Ref. 
HD]). The best agreement is then reached if we addition¬ 
ally assume a dc voltage of U ^ ±5 mV. This is likely the 
voltage offset of the used high-precision current amplifier 
(input) which slowly drifts in time and is hard to control 
due to its dependence on the ambient temperature. Near 
the symmetry points of = 0, the symmetry of the 

overall current is markedly reduced by the dc-current GV 
in Eq. point (iv) above HD]. 

In Fig. H: a) we present measured versus predicted 7(0)- 
curves at Ueff — 0 for two different modulation ampli¬ 
tudes. For accuracy we have introduced the effective 
voltage = U + Vqs corrected by the voltage offset, 
U)ff, caused by the current amplifier. At 0 = tt, Fp and 
Fr oscillate in anti-phase and the QD states Mn are static 
(for all k and independent of A). The modulation of 
Mn(A) grows with |0 — 7r| and with it the number of QD 
levels contributing to current, i. e. fulfilling Mn(A) = Ml,r 
twice during each modulation period. Consequently, be¬ 
tween 0 < 0 < TT (equally for tt < 0 < 27r) |/(0)| in¬ 


creases whenever another QD level starts to contribute 
to For the smaller amplitude at most two levels 

with n = —1,0, for the larger amplitude at most eight 
levels with —4 < n < 3, reach Mn = Ml,r during each 
period and, hence, contribute to 7. As 0 approaches 0 
or 27t the gate modulation becomes symmetric restoring 
TRS, hence 7 drops to zero HD]. The differences between 
theory and experiments in Fig. [^a) are not fully under¬ 
stood, but could be related to the contribution of excited 
states, which are not accounted for in our model. 

As is increased towards |eV^ff| > Eq the dc contri¬ 
bution of the current [first term of Eq. (§] rapidly 

gains weight yielding a transition to axial mirror sym¬ 
metry reflecting the phase independent current direction 
determined by the sign of Ueff. This is clearly visible in 
Eig. I^b) where we include data for the larger amplitude 
at finite V^ff ±0.2 mV. The approximately quadratic 
current increase (decrease) away from 0 = tt (or equiva¬ 
lently from 0 = 0, 27r) is related to the details of Fp .rW- 

The symbols in Fig. [^c) are the measured X-Rs char- 
acteristics at the current maxima near black arrows in 
Fig. lib). Here |Ueff| < E^l2e so that only a single QD 
state contributes to 7 at any instance of time. In this 
regime we find approximately 7 oc Straight lines are 
the corresponding theory data which intersect as a com¬ 
mon point at I4ff 7 ^ 0. The theory data are composed 
of three contributions shown as bar charts: G in blue, 
Qcyciej black and in red. The fair agreement be¬ 
tween theory and experiment confirms the theoretically 
assumed adiabaticity for our experiments. 

In summary, we have realized a Lissajous rocking 
ratchet by modulating a single QD with two phase-locked 
voltage signals at commensurable frequencies. The di- 
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rectional motion is a consequence of breaking the TRS 
which, as our analysis revealed, is restored only for cer¬ 
tain phases between the signals. To achieve this, it is 
not sufficient to modulate the QD-lead couplings but 
additional modulation of the QD levels via capacitive 
cross-couplings is essential. While our experiments used 
a semiconducting QD, similar rocking ratchets could be 
realized in different systems, such as superconducting cir¬ 
cuits, nano-electromechanical systems or molecular elec¬ 
tronics. Lissajous ratchets encode the relative phase, fre¬ 
quency and amplitude information of two radio frequency 
signals into a time averaged dc signal, in spirit similar to 
a Lock-In amplifier. The ability to compare rf-signals 
on the chip promises a refined level of control desired 
for applications related to on-chip spectrum analyzing or 
quantum information processing. 

We thank Klaus Pierz for supplying the wafer ma¬ 
terial. We are grateful for financial support from the 
DFG via LU 819/4-1, SFB-631 and the Cluster of Ex¬ 
cellence “Nanosystems Initiative Munich (NIM)” and by 
the Spanish Ministry of Economy and Competitiveness 
via grant No. MAT2011-24331. S.L. acknowledges sup¬ 
port via a Heisenberg fellowship of the DEC. 


Additional experiments 
Experimental setup 

The CaAs / AlCaAs wafer containing the QD sample 
was mounted on a radio frequency (rf) sample holder as 
shown in Eig.[^ such that the wafer was in direct contact 



FIG. 4. Radio frequency sample holder with sample and a 
scanning electron microscope image of its surface (same image 
as in Fig. [^a) of the main article). 

with a gold plated copper surface thermally connected to 
the mixing chamber of a dilution refrigerator at a base 
temperature of 50 mK. The control gates used to define 
the QD barriers (yellow gates in the detailed sample view 
in Fig. 1^ the gray gate is grounded) are connected to 
stainless steel rf coax cables. Source and drain contacts 
are connected via standard constantan wires and low pass 
filtered at room temperature. All cables are heat sinked 
at several points in the cryostat. A simplified sketch of 
the device circuit highlighting the capacitive control is 
presented in Fig. 

Calibration measurements 

Our numerical calculations rely on accurate measure¬ 
ments of the QD characteristics including the relevant 
energy scales, capacitive coupling constants, tunnel rates 
and realistic line broadening. In the following we present 
our calibration measurements. We begin by showing in 
Fig.i a typical single dot stability diagram depicting the 
current / through the QD as a function of the voltage 
applied to the left gate W and the source-drain voltage 
V. Note that varying the right instead of the left gate 
would result in a comparable plot. The white dots along 
the V = 0 horizontal indicate the positions of current 
maxima in Fig. of the main article. As expected, they 
are located at the crossings of the Coulomb diamonds 
at E = 0. From the dimensions of the Coulomb dia¬ 
monds (horizontal and vertical arrows) we estimate the 
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FIG. 5. Sketch of the QD circuit. Tunnel rates, rL,R, and 
dot levels fin (horizontal lines) are capacitively controlled by 
the gate voltages, Vl and Vr. The leads contain a degenerate 
electron system with chemical potentials, /iL,R- The current 
is measured at the grounded right lead, while the voltage 
V = (/iR — is applied on the left lead. For 1/ > 0 

electrons tunnel via QD levels with /xr > /in > /xl from the 
right to the left lead (blue arrows). 



FIG. 6. Galibration of charging energies and gate voltages 
to energy conversion. Gate voltages are not modulated. Plot¬ 
ted is the current through the QD measured as a function 
of gate voltage Vl = Fl + 202 mV and bias Vr at constant 
Fr = Fr + 91 mV = —2.5 mV. The numbers indicate our in¬ 
dex parameter n. White dots mark the corresponding current 
maxima positions in Fig. of the main article. Red lines il¬ 
lustrate the Goulomb blockade diamonds. 


charging energies Ec{n) = /Cqj:,{n) and the capaci¬ 

tive coupling = eCL,R(n)/CQD(n), where Cqj:>{n) is 
the eigen capacitance of the QD [at the working point 
Cl — C^ 50 aF and Cql) — 130 aF]. We find an ap¬ 
proximate linear dependence Ec{n) = Eq + n5E with 
5E = —0.04 meV and at the working point Eq = 1.27 
meV (solid lines in Fig. and an approximately con¬ 


stant leverage factor = 0.38e. Thus, the gate volt- 
ages^shTt the chemical potentials of the resonances as 
Ain(rL, Ui) = +eo(VL, Vr) with = /U„(0,0) and the 
onsite energy 


eo{VL,VR) = -K^{VL + VR). (3) 

The above relations can be reinterpreted as /x^ = nE^ — 
n{n-\-l)SE/2-\-jj.^{n = 0), where /x^(n = 0) = 0.41 meV is 
the difference between the chemical potential of grounded 
leads ajid that of the lowest unoccupied dot-level at 
= 0. The influence of the gate voltages on 
the dot-lead tunnel rates can be estimated from the dc 
current shown in Fig. [^a) of the main text. In Fig. 
we plot these data along the bisecting line Fl = of 



FIG. 7. Galibration of tunnel couplings. Gurrent through 
the QD for F = 100/xV, data from Fig. along its bisecting 
line where Fl = Fr and Fl — Fr (black dots). Also shown 
are two theory curves, all using Kr = 0.03e, F = 0.13 meV 
in Eq. charging energies from Table and assuming 

Fl = Fr: eFL,R/h according to Eq. Q current calculated 
with standard scattering theory with finite tunnel barriers as 
only broadening mechanism (red solid line), see Eq. and 
below; calculated current as above but additionally convo¬ 
luted with a Gaussian inhomogeneous distribution of constant 
width mimicing slow charge noise (blue line). 

Fig. a). The Coulomb blockade oscillations display cur¬ 
rent maxima, but their smoothed average exhibits an ex- 
ponenDal d^endence on gate voltages below the working 
point Fl = f^R = 0 devolving into a linear increase above. 
To capture this averaged behavior with only two free pa¬ 
rameters, we assume the dependence as 


TaiVa) 


_ l^rVg _ 

1 - ex.p{-KrVa/T) ’ 


( 4 ) 


where a = L,R, while F = F^ = = 0.13 meV is the 

tunnel coupling at the working point and Kr — 0.03e is 
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TABLE L Experimentally determined parameters of our 
QD; from left to right: index number n, chemical poten¬ 
tials at the working point El = ^ = 0, charging energies 
Ec{n) = l-in— tunnel couplings for the case of symmet¬ 

ric coupling El = Er and El = Er such that fin = 0, i.e., 
Tl.r = ra(/i“/2K^), see Eqs. § and 0. 


n 

(meV) 

Ec(n) (meV) 

El,r (meV) at /in = 0 

-3 

-3.85 

1.39 

0.39 

-2 

-2.25 

1.35 

0.59 

-1 

-0.90 

1.31 

0.82 

0 

0.41 

1.27 

1.09 

1 

1.68 

1.23 

1.40 

2 

2.91 

1.19 

1.74 

3 

4.11 

1.15 

2.13 

4 

5.31 

1.11 

2.53 



the slope for large R, Lr,. The red solid line in Fig. 
represents the prediction from Eq. © below based on 
scattering theory including the finite couplings El(E) = 
rR(^) as only broadening mechanism [discussed below 
Eq. Convoluting this curve with an additional Gaus¬ 
sian broadening of constant width yields a better fit to 
the measured data (blue line). The Gaussian expresses 
an inhomogeneous broadening caused by slow statistical 
fluctuations of the QD levels (charge noise). In our case 
this inhomogeneous broadening is the dominant broad¬ 
ening mechanism. The non-perfect fit is likely a con¬ 
sequence of the increasing relevance of dissipation as the 
barriers open up. Note that the values of El(E) = Er(E) 
are predetermined by the integral of the current peaks 
which is independent of additional inhomogeneous broad¬ 
ening. In Table |T] we summarize the functional depen¬ 
dences along the bisecting line in Fig. [^a) of the main 
article according to the equations above. 


Radio frequency calibration — proof of adiabaticity 

Next, we discuss the calibration of the amplitude of 
the radio frequency gate voltage modulation. Figure 
demonstrates an apparent problem of our non-perfect 
rf setup, namely strong cable resonances, i. e. standing 
waves caused by reflections of the rf-signal at the sample 
holder and meters away at the rf generator. In the adi¬ 
abatic regime we expect to find a pump current |/| oc /, 
if E = 0. Instead we observe strong current oscilla¬ 
tions which indicate an oscillating modulation amplitude 
even though the rf signal strength is fixed, the main rea¬ 
son being cable resonances. For our ratchet experiments 
we restricted ourselves to specific modulation frequen¬ 
cies, namely / = 25, 50,100,150, 200 MHz. In Fig. we 
demonstrate exemplarily for one of the gates the cali¬ 
bration of the modulation amplitude at these frequen¬ 
cies. Plotted is the differential conductance of the QD as 
function of Vr and its modulation amplitude (where 


FIG. 8. Maximum current /max for phase difference cj) from 
0 to 27r through the unbiased QD plotted as a function of 
modulation frequency at A; = 1 and modulation amplitude 
A = 3/;g. 


7.5 

> 

-7.5 


/= 25 MHz 


50 MHz 


-7.5 



8.2 


- 8.2 



0 Lf(mV) 


Fr (mV) 400 


FIG. 9. Radio frequency calibrations of the right hand side 
gate for / = 25, 50,100, 200 MHz in panels a-d. Plotted is the 
differential conductance g — dl/dV as a function of applied 
modulation voltage Er^ and mean gate voltage Vr (at con¬ 
stant El = —90 mV). The flocks of parallel black lines with 
mutually opposite slopes are guide for the eyes indicating cur¬ 
rent resonances which depart proportional to the increasing 
amplitude. 


ER(t) = Er + ER^cos 27 r/t) while the left gate voltage 
is kept constant. Each Coulomb resonance splits in two 
resonances located at the turn around points of ER(t). 
These splittings are proportional to and are indicated 
by black lines in Fig.^ Their slopes cxR{f) = Vr/V^ can 
be used to calibrate the actual rf modulation amplitude 
























TABLE IL Radio frequency calibration factors aL,R at vari¬ 
ous driving frequencies. 


/ in MHz 

1/ol 

l/OiR 

25 

25.2 

24.9 

50 

23.2 

23.1 

100 

27.5 

27.2 

150 

26.3 

26.1 

200 

182.5 

178.6 


according to: 

^L,R = ^M<^L,R(/)^LfR/e (5) 

The result of the calibration (of both gates) is summa¬ 
rized in Table [nj To account for the frequency dependent 
calibration in our measurements we correct the ampli¬ 
tudes at each frequency according to Table [ll| and Eq. 

In Fig.j^a) we present /((/)) at four frequencies for /c = 1 
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EIG. 10. (a) Current / as a function of phase difference 0 

for /c = 1, and A = 3Eq. The current has been corrected 
for the frequency dependent transfer functions of the cables 
by applying the calibration factors ol,r(/) listed in Table [n| 
according to Eq. (b) Same data as in panel a but scaled 
by/(/)^ ^^»/(/). 

and identical modulation amplitudes A = 3Eq after cal¬ 
ibration. The same data scaled by /(/) ^ ioomhz 
to the current expected at / = 100 MHz are approxi¬ 
mately frequency independent [Fig. [TQ|[b)]. It demon¬ 
strates not only the validity of our procedure to calibrate 
the rf amplitudes but also corroborates our assumption of 
adiabatic charge transport. Small deviations, especially 
those near 0 = 0 are probably related to a local disorder 
potential which compromises the spacial symmetry. The 
importance of a correct calibration becomes evident in 
Fig. [TT| plotting /max in direct comparison to the uncali- 
brated current taken from Fig.[^ Straight lines are guide 
for the eyes indicating the adiabatic transport regime. 



FIG. 11. Frequency dependence of current before and af¬ 
ter calibration of the modulation amplitude. The dashed 
line presents uncalibrated data identical to those in Fig. 
Golored symbols indicate current values of specific maxima 
and minima (constant phase) marked in Fig. 10 a) by arrows; 
same as in the inset of Fig. ^b) of the main paper. The 
straight lines are a guide for the eyes and verify the expected 
adiabatic regime. 


Theoretical approach 


We derive a scattering formalism for the dc current 
through a periodically driven conductor in the low- 
frequency limit, which generalizes the scattering ap¬ 
proach for parametric pumping [TSl [16] to the presence 
of a finite dc voltage and a tiny ac bias. The derivation 
starts from Floquet scattering theory [laiTTHis] which 
we expand for a small driving frequency / = Q/27r. De¬ 
spite that we will finally apply the approach to a two- 
terminal device, we formulate the general theory for an 
arbitrary number of leads. 

A cornerstone of Floquet scattering theory is an ex¬ 
pression for the probability Tc^/ 3 (t, e) that an electron with 
initial energy e from lead 0 is scattered to lead a at time 
t. Then the time-dependent net current in lead a can be 
expressed as muis] 


/< 




/(e - - FI Tapit, e)/(e “ /i/3(^)) 


(6) 

where h = 27Th denotes Planck’s constant while f{x) = 
[exp(x//c^T) + 1]“^ is the Fermi function for thermal 
energy ksT. In our experiment, the modulation of the 
gate voltages induces a tiny ac bias caused by a tiny ca¬ 
pacitance between the modulated gates and the electrical 
leads of the QD. To capture this, we allow for periodically 
time-dependent chemical potentials fia {t) = {t) • 

Our goal is to derive a convenient expression for the time- 
average of Ia{t) in the adiabatic limit. 
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We start by Taylor expansion of the Fermi functions in 
the ac shift Wa{t). It yields the time-dependent current 
Ga^{t)wi 3 {t) with the time-periodic conductance matrix 

Gapit) = Ga0(v{t)) = ^ [Sa/3 “ fl)] . (7) 

Despite that the ac bias Wjg(t) vanishes on average, its 
impact becomes rectified and contributes 

dt 

= I ^GMmMt) ( 8 ) 

to the dc measured current, where T is the modulation 
period. 

Having considered the ac bias, we now come to the 
effects of a dc voltage and the modulation of tunnel bar¬ 
riers and the QD levels. We employ ideas put forward 
in Refs. HHiEniEi]. For the now constant chemical po¬ 
tentials, the time-average of the transmission TQ,^(t,e) 
provides which is the probability that an elec¬ 

tron with initial energy e coming from lead (3 is scattered 
to lead a after the absorption of k energy quanta, where 
negative k corresponds to emission. Thus, the final en¬ 
ergy of the electron is e + khCl. Then the net average 
current flowing from lead (3 to the conductor assumes 
the intuitive form [laiMin] 


I. = 



/3,k 



( 9 ) 


This expression can be derived from a time-dependent 
wire-lead model which relates the transmission to the 
scattering matrix and the Green’s function of the driven 
conductor [la na Eg. From the unitarity of the scat¬ 
tering matrix follows charge conservation, which means 
that an electron with initial energy e will be scattered 
with probability unity to any sideband k in any lead a, 
so that the sum over all these channels must fulfill the 
relation 


a,k 

With the same arguments but summing over all initial 
states that end up with energy e in lead a we obtain a 
further sum rule: 


Eq. (© as prefactor of /(e —/io,) into the current formula 
^ and obtain 

4 = E ^ - fie - Ai/3)] 



Notice that in contrast to Ref. [18], we have expanded in 
O for the transmission rather than for the Fermi function, 
so that the result holds whenever the transmission as a 
function of e is smooth on the scale hCl. In particular, 
it is valid also for low temperature, i.e. beyond the high- 
temperature limit of ksT hfl. 

Next we make use of the time-periodicity of the driving 
which has the consequence that the propagator U— 
depends not only on the difference r, but explicitly on 
both the initial time t — r and the final time t. After 
Fourier transformation with respect to r, one obtains 
both the propagator and the scattering matrix for an 
electron with initial energy e [12]. The result is periodic 
in time and can be written as 


Sit,e) = Sit + T,e) = E (13) 


Then from the usual relation between transmission am¬ 
plitudes and probabilities, e)p, and 

taking the time-average of the current follows (e) = 
[Notice that generally the transmission to the 

kth sideband, (e), is different from the kth Fourier co¬ 
efficient of Tap{t, e)]. Next we employ Parseval’s theorem 
to write the /^-summation in Eq. (12) as time-integration 


over one driving period. Moreover, it is convenient to de¬ 
fine voltages Va as deviations from an average chemical 
potential, /io, = /i + eW, and evaluate the Eermi func¬ 
tions for zero temperature and small voltage. Then we 
end up with the time-averaged current 




Qcycle 


(14) 


with the time averaged conductance 

Gap ~ Y ^"/ 5 (^(^)) 


(15) 


J2T3pie-khi}) = l. ( 11 ) 

(3,k 


Adiabatic limit 

We have performed experiments in the adiabatic limit 
corresponding to an expansion of la to first order in Cl. 
Eor the Taylor expansion we insert the left-hand side of 


and the current resulting from the rectified ac bias given 
in Eq. (§. The last term in Eq. ( [T^ reflects the charge 
parametrically pumped through contact a during one 
driving period. 




The shape of the driving can be expressed by a closed 
curve C in parameter space, v{t) = v{t + T). Then, if C 
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is traversed adiabatically slowly, the scattering matrix 
depends merely parametrically on time, i. e. 5'(t,e) = 
5'(iT(t), e), where the latter is the static result for the 
instantaneous value of v at time t. This allows one to 


transform the time integral in Eq. (16) to a closed line 
integral to obtain 




Applying Stokes’ theorem, one can transfer this expres¬ 
sion into a surface integral to obtain Brouwer’s formula 
m- For the numerical evaluation in the parameter range 
considered here, however, the present form is more ap¬ 
propriate. 


An important implication of Eq. (14) is the separa¬ 
tion of the average current into a dc contribution stem¬ 
ming from the average conductance and an adiabatically 
pumped charge. The main difference between these two 
quantities is their behavior under time inversion of the 
closed curve C: While the average conductance is invari¬ 
ant, the pumped charge acquires a minus sign. Below, 
we will explore this symmetry property for the Lissajous 
curves applied in the experiment. 


Independent channel approximation 

In our experiment, the conductor is formed by a QD 
connecting two leads, a = L^R. It is driven with an 
amplitude A ranging from zero to the rather large value 
5 Eq , so that the onsite energy of the QD may change by 
several charging energies—we will consider up to 10 ex¬ 
cess electrons. Moreover, the influence of the gate volt¬ 
age applied to the tunnel barriers turns out to be cru¬ 
cial for the observed pumping, so that a master equation 
approach based on lowest-order perturbation theory in 
the dot-lead coupling is not appropriate. As a conse¬ 
quence, already for the static transport problem, a full 
treatment that includes all possible spin and correlation 
effects is practically impossible. We instead assume that 
each charge state of the quantum dot contributes as inde¬ 
pendent transport channel to the current. This requires 
that the separation of the resonances given by the charg¬ 
ing energy is larger than their homogeneous widths such 
that the conductance peaks are well separated. More¬ 
over, the driving must be slow enough that the quantum 
dots can be assumed to be always relaxed to its many- 
particle ground state. Our QD model is parametrized by 
the energies fin necessary for adding the next electron 
to the QD. Thus, the position of the corresponding con¬ 
ductance peak is directly related to Assuming that 
all these peaks are well separated, we can treat the res¬ 
onances as independent levels that can be occupied by a 
spinless electron. The corresponding scattering matrices 


read 


Sn{e) = 1 - 


_ i _ ( Tl VriTV\ 

e - Mn + «(rL + rR)/2 Va/TlTr Tr )' 

(18) 

where Tl and Tr are the dot-lead tunnel rates. 

We drive our QD by applying time-dependent gate 
voltages to the tunnel barriers. Thus these gate volt¬ 
ages form our parameter space and are the components 
of the vector = [^(t), ER(t)], see Eq. (1) of the 

main article. They entail a time-dependence on the on¬ 
site energy eo and thus on as well as on the dot-lead 
tunnel rates Tl and Tr. For the onsite energy, we as¬ 
sume the linear relation The height of each tunnel 
barrier is also shifted approximately linearly by the re¬ 
spective gate voltage. For ideal tunnel barriers according 
to the WKB formula, one would expect that the dot-lead 
tunnel rates depend exponentially on the gate voltages. 
Our dc measurements, however, indicate that this over¬ 
estimates the rates for large positive El,r- Therefore, we 
assume for rQ,(fT) the dependence given by Eq. with 
T = Tl = = 0.07 meV and Kr — 0.03. In order to 

evaluate the current formula ( p^ , we insert the scatter¬ 
ing matrix Sn into Eqs. (15) and (p^Tl) and subsequently 


sum the contribution of the resonances n. 

Even though this approach already allows us to obtain 
the dc current (all our numerical results are computed in 
this way), it is instructive to investigate the contribution 
of a single resonance to both the average conductance 
and the charge per cycle. For this purpose, we assume 
that at a resonance, the dot-lead rates are weakly time- 
dependent with average values T/^^r, while /in(^) — {t — 
to)fin^ Then the straightforward evaluation of Eq. ( p!^ 
with the off-diagonal matrix element of Sn yields 


- 


TlEr n 
^ Tl + Tr I/in 


(19) 


Notice that /in is proportional to Q and, thus, G'^ is 
frequency independent. 

The pumped charge is essentially determined by the 
scattering phase in the prefactor of S'n- Assuming again 
that the Tare practically constant at the resonance, 
we find [21] 




efL 


bn + Tr 


sign(/in). 


( 20 ) 


This means that for Tl Tr, an electron enters from or 
leaves to the left lead depending on the sign of /in- In the 
opposite limit, the right lead is relevant and QJf^ is much 
smaller than the elementary charge. This also emphasizes 
the role of the parameter dependence of the tunnel rates: 


In the derivation of Eq. (20), we assumed that the Tl 


stays constant when a level crosses the chemical potential 
of a lead at time t = to. However, in a complete pump 
cycle, the level will cross the chemical potential as many 
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FIG. 12. (a) Current in the absence of the modulation for bias voltage V = 200mV. (b) Component ai(v) of the vector field 

(21) appearing in the symmetrized line integral representation of the charge pumped per cycle. Positive values (blues stripes, 
marked with +) emerge when a QD level equals the chemical potential of the left lead and correspond to electron pumping from 
the left lead to the QD. Accordingly, the curvature is negative (red, —) when a QD level agrees with the chemical potential 
of the right lead, while the mismatch stems from the voltage bias. In the relevant regime, the second component, a 2 (v), looks 
very similar, (c) Corresponding curvature u; = (daxjdx^ — da^jdvx) that appears in the symmetrized Brouwer’s formula. The 
lines are splitted due to the bias. 
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FIG. 13. Theoretical prediction of the current for k = 1 shown in Fig. [^b) of the main article, /l = /r = 200 MHz. The dc 
tunnel current GV (a), the pumped current (b), and the ac current (c) sum up to the total current plotted in (d). 




FIG. 14. Theoretical prediction of the current for k — 2 shown in Fig. |^a) of the main article, /l = 50MHz /r = lOOMHz. 
The dc tunnel current GV (a), the pumped current (b), and the ac current (c) sum up to the total current plotted in 

(d). 


times from above as it crosses from below. Thus if Fl,r 
were constant, the net charge pumped from the left lead 
to the dot in the whole cycle would vanish. In turn, we 
will observe significant pumping from the left lead to the 
right lead if Fp ^ Fr when the dot level is lowered, while 
Fp Fr when the level is raised. 

For a two-terminal device, the pumped currents enter¬ 
ing the dot from the left and from the right lead compen¬ 


sate each other so that + = 0. Nevertheless, 

the vector fields appearing in the integral in Eq. 0 may 
differ by more than their sign, because the integrands are 
relevant only up to a gradient field which vanishes under 
the closed line integral. Therefore it is more convenient 
to discuss the integrand of the symmetrized expression 
for the pumped charge: = 
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dv • a{v) with the vector field 


a = 


is 


{(S.gradSD^j^ -(S.gradS*)^^! }. 

( 21 ) 

Figure [T^b) shows the first component of d{v) for the 
parameters of our quantum dot. A significant contri¬ 
bution requires that the resonant level lies close to the 
chemical potential of a lead. For small bias, this is also a 
necessary condition for the emergence of a conductance 
peak, see Fig.j^a). To make direct use of Brouwer’s for¬ 
mula, one would transform via Stokes’ theorem the line 
integral into a surface integral [see remark after Eq. 0] 
and would end up with an integration over the curvature 
shown in Fig. fl^c). 


Lissajous pumping and symmetries 

The gate voltages applied in the experiment follow Lis¬ 
sajous curves in parameter space defined as 


••”=( 11 !)- 


= A 


cos(Ot — 0) 
cos{kQt) 


( 22 ) 


where henceforth k is referred to as order and 0 is the 
phase shift between the two components. Driving the two 
barriers with different but commensurate frequencies is a 
particular feature of our experiment. The amplitude A is 
equal at both barriers and is typically several Eq in our 
experiments. A central question is that of the symmetry 
properties of the dc conductance and the pumped charge 
as a function of the modulation parameters. 

We start our symmetry considerations by noticing that 
the Lissajous figures obey Vk,cf>{t) = f4,(/)+27r//c(^+27r//cr^), 
namely that a phase shift of 27r/k is canceled by a time 
shift of 2711kft. This implies that the dc current patterns, 
being averaged over time, as a function of 0 possess a k- 
fold symmetry in <p. 

A further symmetry property of the three contribu¬ 
tions to the current in Eq. (14) follows from the behavior 


of Ga/3: and under time reversal. It plays 

a crucial role for adiabatic pumping m and ratchet ef¬ 
fects. Below, we find 2k symmetry points at which Ga^ 
is symmetric, while and are anti-symmetric 

under time reversal. Eor a derivation, we determine for 
each Lissajous curve Vk,cf) a time-reversed partner with 
phase (p' which must fulfill the condition 




(23) 


for all times t. The inversion t ^ — t is thereby accom¬ 
panied by a time offset to, which is permitted by the 
time periodicity of the integrands in Eqs. ([^, ([T^, and 
( [T7[ ). Inserting Eq. ( [^ into Eq. (23) yields ftto = 2711 jk 
and p' = —p + 27\(ijk + i') where t = 0,1,..., /c — 1 


and E = 0,1. The special phases pggi = 7r(-t//c F E) 
fulfill p = p' (with to given above), so that the origi¬ 
nal and the time-reversed Lissajous curves [defined by 
Eq. (22)] he on top of each other, while they evolve in 
opposite direction in time. Thus, the phases pii> de¬ 
fine 2k in-equivalent points with TRS. An interesting 
observaticm is that there, the Lissajous curve takes the 
form Vr(\A) = ±cos[/carccos(±kL)]- This defines the 
kth Chebyshev polynomial [23] which represents a de¬ 
generate loop that does not enclose a finite area. As a 
consequence, and F^^^ (with odd symmetry under 

time-reversal) vanish. Hence, at the symmetry points the 
current takes the value GV and, in particular, vanishes 
for E = 0. 


Substituting in Eq. (|22|) the phase p by the deviation 

, one can see 
—Ap. Thus, 


from the symmetry point, Ap = p — 
that time reversal corresponds to Ap 
the behavior of the integrals in Eqs. (|^, (15), and ( [l7| ) 
under time reversal allows us to draw conclusions about 
the symmetries of G^^ and as a function of Ap. 

To be specific: Since the average conductance Eq. (15) is 


invariant under time reversal, it must be an even func¬ 
tion of Ap. By contrast, the charge pumped per cycle, 
Eq. ( [T7| ), changes its sign upon time reversal. Conse¬ 
quently, it must be an odd function of Ap. In short, 
despite the only k-fo\d symmetry of the Lissajous curves 


Eq. (22), we find 2k in-equivalent phases at which Ga^ 
is symmetric, while is anti-symmetric. 

Our experimental results indicate that the rf mod¬ 
ulation of gate voltages induces a tiny ac bias. Its 
origin is a tiny capacitive coupling between the mod¬ 
ulated gates and the 2D leads of the QD. Therefore, 
we model it as modulation of the chemical potentials 
/^a(^) = + ewa{t) with a phase shift of —7r/2 com¬ 

pared to the modulation of the gate voltages [described 
in Eq. (1) of the main article]: 




— p) 
sm{kflt) 


(24) 


The capacitive coupling constant ACac turns out the be of 
the order 5 • 10“^ and is expected to decrease with the 
modulation frequency, /^ac 


In comparison to the gate voltages in Eq. (22), the ac 


modulation of the chemical potentials, uJ, contains a sine 
instead of a cosine. In a symmetry analysis along the 
lines above, this finally leads to a minus sign. Therefore 
the integrad in Eq. <© is an odd function of t (besides 
a time shift by to), so that the rectified current changes 
its sign under time reversal. Consequently, F^^^ is anti¬ 
symmetric in Ap. 

In Eigs. [T^and[T^ we show for /c = I, 2 the three con¬ 
tribution to the current in Eq. <§ of the main text, i.e., 
the dc current GE, the pump current and the 

contribution of the rectified ac bias, F^^^. This visualizes 
the symmetry of the dc current and the anti-symmetry 
of the second and the third contribution. Notably, for 
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the case of a tiny dc voltage ofV = l/SV considered here, 
the overall behavior [panel (d)] is dominated by the anti¬ 
symmetric contributions. 


Additional data 

In Fig. [^we compare the measured Lissajous ratchet 
data for /c = 1, 2, 3 also discussed in the main article but 
all in spherical coordinates and add another data set for 
k = A. The data contain the predicted k -1- 1-fold symme¬ 
try but also show a complex detailed structure caused by 
the combination of the three contributions in Eq. of 
the main article, namely the time averaged conductance 
GV, the ratchet current and It is possible 

to reduce the information in similar measurements, such 
that the k + 1-fold symmetry is more evident, by increas¬ 
ing the dot-lead tunnel couplings. We do this by shifting 
the working point to the mixed valance regime where the 
Coulomb current peaks overlap strongly (W = 12.5mV 
and Vr = 12.5mV, see Fig. Ba) of the main article). The 
result is a strong broadening which thoroughly washes 
out most fine structure of the current. In Fig. we 
present such data up to /c = 3. These measurement are 
a clear demonstration of the Lissajous ratchet effect as 
they emphasize the symmetry properties of 


Applications and alternative realizations 

In our Lissajous rocking ratchet a combination of rf 
excitations breaks the time reversal symmetry and yields 
a dc current. It represents a new method to access radio 
frequency information by mutually comparing the fre¬ 
quencies, amplitudes and phases of rf signals. As such, 
Lissajous rocking ratchets could find applications as rf 
comparator or serve as detector or filter components. 
Our quantum dot implementation is interesting for quan¬ 
tum information processing applications, as it allows one 
to compare and process small rf-signal on-chip. In a sim¬ 
plified picture a Lissajous ratchet resembles an on-chip 
rf Lock-In amplifier: to illustrate this idea imagine that 
we apply a clean rf reference to one gate and a noisy sig¬ 
nal to the second gate. The resulting dc current strongly 
depends on the relative phase, amplitude and frequency 
of the two signals. By measuring a time averaged dc 
current, our device filters the time periodic carrier signal 
while averaging out uncorrelated noise, just as a Lock-In 
amplifier does. However, the Lissajous ratchet allows a 
detailed comparison of mutual rf-signals going well be¬ 
yond the scope of a Lock-In amplifier. 

Lissajous rocking ratchets are not limited to QD cir¬ 
cuits but could be realized in a variety of systems, for 
instance in macroscopic electronics or mechanics or in 
superconducting circuits. In the following we propose a 
possible realization in a nanoelectromechanical systems. 


In Fig. 17 we extend the already realized “nano bell” [24] 
and illustrate what we mean by “nanomechanical New¬ 
ton cradle”. The device is composed of two cantilevers 
and a fixed metallic contact (R). The central cantilever 
is patterned with an isolated metallic island while the 
left one contains an island but electrically connected to 
a second lead (L). The islands are positioned such, that 
the central island can touch both the left and right leads 
whenever both cantilevers are vibrating with sufficiently 
high amplitude. A voltage can be applied between L and 
R and current can be measured. We assume that both 
cantilevers can be externally driven (at predefined rela¬ 
tive phases) at their mutual eigen frequencies and that 
the latter are tunable, e. g. , via capacitive coupling to 
additional gates or by using dielectric forces [25] . We also 
assume that the central island is small enough to allow 
for single electron transport based on Coulomb blockade, 
namely that we can treat it as a QD [24] . The electronic 
levels of the islands are then modulated by the same ca¬ 
pacitive coupling that is used to mechanically drive the 
cantilevers. In such a device the coupling between the 
islands and that between contact R and the center island 
are strongly time dependent, zero most of the time and 
strong during touching. This system is highly tunable 
and very close in spirit to our QD circuit: importantly, 
the time dependence of the couplings and island levels 
are linked, a central precondition to define a Lissajous 
ratchet. 
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FIG. 16. Measured current / versus modulation amplitude and phase for (a) /c = 1 and /l = /r = 100 MHz, (b) /c = 2 and 
/l = 50 MHz /r = lOOMHz, (c) /c = 3 and /l = 50 MHz /r = 150MHz. The working point is defined by Vl = 12.5mV and 
Vr = 12.5mV near the upper right corner of Fig. |l[a ) of the main article. The enhanced dot-lead tunnel coupling washes out 
most fine structures compared to the data in Fig. ]l5| above. 
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FIG. 17. Nanomechanical Newton cradle discussed in the 
main text. 








